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A quantum algorithm is proposed to solve the Satisfiability problems by the ground-state quantum 
computer. The scale of the energy gap of the ground-state quantum computer is analyzed for the 3- 
bit Exact Cover problem. The time cost of this algorithm on the general SAT problems is discussed. 
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I. INTRODUCTION 

A quantum computer has been expected to outper- 
form its classical counterpart in some computation prob- 
lems. For example, the well-known Shor's factoring 
algorithm^] and Grover's algorithm^ accelerate expo- 
nentially and quadratically compared with the classical 
algorithms, respectively. It is a challenge to find whether 
a quantum computer outperforms on other classically in- 
tractable problems 0, 0| , which cannot be solved classi- 
cally in polynomial time of N, the number of the input 
bits. 

Especially interesting are the NP-complete 
problems 0, which include thousands of problems, 
such as the Traveling Salesman problemjfj and the sat- 
isfiability (SAT) problems. All NP-complete problems 
can be transformed into each other in polynomial steps. 
If one of the NP-complete problems can be solved in 
polynomial time by an algorithm even in the worst 
case, then all NP-complete problems can be solved in 
polynomial time. However, it is widely believed that 
such a classical algorithm doesn't exist. In this paper we 
will discuss quantum algorithm for solving SAT prob- 
lems. A K-SAT problem deals with N binary variables 
submitted to M clauses with each clause Cj involving 
K bits, and the task is to find iV-bit states satisfying 
all clauses. When K > 2, K-SAT is NP-Complete, and 
some instances become classically intractable when the 
parameter a = M/N, as M, N — > oo, approaches the 
threshold a c (K ) 0, S 0, G3 . 

Due to the properties of quantum mechanics, it's hard 
to design quantum algorithms directly from intuition. In 
the present paper, we will study the properties of the 
ground-state quantum computer(GSQC), and show that 
the special property of the GSQC naturally leads to al- 
gorithm for solving SAT problems. Although we cannot 
determine whether or not this algorithm solves the NP- 
complete problems in polynomial time, we try to shed 
light on the complexity of the NP-complete problems. 

In the following sections, at first we introduce the idea 
of the ground-state quantum computer |lllll2Ul3| and its 
energy gap analysis 14] , then demonstrate the particular 
property of the GSQC, which provides a direct approach 
to solving SAT problems, and finally an example, an al- 
gorithm for solving the 3-bit Exact Cover problem, is 
given. 



II. GROUND-STATE QUANTUM COMPUTER 
AND ITS ENERGY GAP 



A standard computer is characterized by a time- 
dependent state \4>(ti)) = Ui\ip(ti-i)), where U denotes 
the instance of the i-th step, and Ui represents for a uni- 
tary transformation. For a GSQC, the time sequence is 
mimicked by the spatial distribution of its ground-state 
wavefunction |"0o)- As proposed by Mizel et.al.pH- the 
time evolution of a qubit may be represented by a col- 
umn of quantum dots with multiple rows, and each row 
contains a pair of quantum dots. State |0) or |1) is repre- 
sented by finding the electron in one of the two dots. It 
is important to notice that only one electron exists in a 
qubit. The energy gap, A, between the first excited state 
and the ground state determines the scale of time cost. 



A. Hamiltonians of GSQC 

A GSQC is a circuit of multiple interacting qubits, 
whose ground state is determined by the summa- 
tion of the single qubit unitary transformation Hamil- 
tonian h J (Uj), the two-qubit interacting Hamiltonian 
h(CNOT), the boost Hamiltonian h(B,X) and the pro- 
jection Hamiltonian h(\j), A). 

The single qubit unitary transformation Hamiltonian 
has the form 

h j (Uj) = e [C^C^ + C]C 3 - (C]l I', i + /i.e.)] ,(1) 
where e defines the energy scale of all Hamiltonians, 



C 
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~3,0 



i is the electron creation operator 

on row j at position 0, and Uj is a two dimension matrix 
representing the unitary transformation from row j — 1 
to row j. The boost Hamiltonian is 



h 3 {B,\) 



A 



(2) 



which amplifies the wavefunction amplitude by the large 
value number A compared with the previous row at \ipo}- 
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The projection Hamiltonian is 



B. Energy Gap of GSQC 



tf'(| 7 >,A) = 



c j-i,-f C j—h'y + 
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, (3) 



where I7) is the state to be projected to on row j and to 
be amplified by A at \ipo). The interaction between qubit 
a and (3 can be represented by h(CNOT): 

hi^CNOT) 

+4,^,^(1) + ct dA c a>j>1 h^(N). (4) 

where for c\ b , its subscription a represents for qubit a, 
b for the number of row, 7 for the state I7). With only 
hP{Uj) and h^JCNOT), its ground state isp3: 
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All above mentioned Hamiltonians are 
inite, and are the same as those in 
pairwise interaction is considered. 

The input states are determined by the boundary con- 
ditions applied upon the first rows of all qubits, which can 
be Hamiltonian h = E(I + <2j<7i) with cr^ being Pauli 
matrix and J^. af = 1. For example, with h = E(I+a z ), 
\ipo) on the first row is |1); with h° = E(I — <r x ), it is 
(|0) + |1)). If E is large enough, for example, at E > lOe, 
the energy gap will saturate and become independent of 
the magnitude of E [l4j . 

To implement an algorithm, on final row of each qubit 
a boost or a projection Hamiltonian is applied so that 
IV'o) concentrates on the position corresponding to the 
final instance in the standard paradigm, hence measure- 
ment on the GSQC can read out the desired information 
with appreciable probability. With boost Hamiltonian 
or projection Hamiltonian on last rows, the ground-state 
wavefunction amplitude on those rows will be A of that 
on their neighboring rows. 

By observing the expression Eq.|J5J, it's easy to find 
that, for two interacting qubits, the ground-state wave- 
function has the form 1141 
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where each qubit is divided by the interacting Hamilto- 
nian as two parts, and the part with boundary Hamilto- 
nian h° is called as upstream, and the other part is called 
downstream. In this paper, we always use this definition 
when upstream or downstream is mentioned. 



Now we briefly introduce how to find the scale of the 
energy gap of a GSQC. For details, please find in [l4j . 

With multiple interacting qubits, one needs to evaluate 
on each qubit the parameter 1/x, the overall amplitude of 
lowest excited state on top rows of this qubit before meet- 
ing the first interacting Hamiltonian, assuming that on 
the top rows of this qubit the lowest energy excited state 
is orthonormal to |^q) while states on all other qubits 
remain the same as the corresponding ground state with 
only magnitude changed. The energy gap|l4j is given by 
the minimum parameter 1/x as 



A cx e(l/x) 
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(7) 



The rule of estimating 1/x is as following [l4j : With 
each qubit ended with either a projection or a boost 
Hamiltonian containing the same (for simplicity) ampli- 
fying factor A ^> 1, when estimating 1/x for a qubit, say 
qubit A, (i) at first x is set to 1; (ii) the boost Hamilto- 
nian, not the projection Hamiltonian, on qubit A itself 
increases x by multiplication of A; (Hi) if qubit A directly 
interacts with another qubit, say qubit B, by Hamilto- 
nian tiAB, then we determine, excluding qubit A, on the 
qubit B the ground-state wavefunction amplitude ratio 
of the upstream part (with respect to Hab) over its final 
row, — , contributions to — are found one by one ac- 

' XB ' XB J 

cording to Eq. @ : if the upstream part of qubit B doesn't 
coexist with the states on final rows of any one qubit, ex- 
cept for qubit A, then xb should be multiplied by a A; 
(iv) finally, the value of 1/x on qubit A should be mul- 
tiplied by or Ilj— if more than one qubit directly 
interact with qubit A. 

According to the above rule, the energy gap A of single 
qubit with length n and ended with boost Hamiltonian 
h(B,X) scales as e/A 2 as A > n; when ended with pro- 
jection Hamiltonian /i(|7), A), A is independent of A. For 
two n-row qubits interacting by h(CNOT), A cx e/A 4 
as A 3> n if both qubits ended with h(B, A) or one with 
h(B,X) and the other with h(\j),X). Numerical calcula- 
tions confirm these results. The Fig. (lb) and Fig. (2) in 
|l4j are two examples on how to apply the above rule on 
complicated circuits. 

Complicated GSQC circuit may have exponentially 
small energy gap, like the circuit in Fig. (lb) of 01 > an d 
assembling the GSQC circuit directly following the al- 
gorithm for the standard paradigm, such as quantum 
Fourier transform, leads to exponentially small energy 
gap. In order to avoid such small gap, the teleportation 
boxes are introduced on each qubit between two control 
Hamiltonians Q. Fig.QJ shows how the CNOT interact- 
ing qubits is modified by inserting teleportation boxes on 
each qubit 's upstream and downstream part. The tele- 
portation boxes make all qubits short (the longest qubit 
has length 8), on the other hand, for arbitrary GSQC cir- 
cuit they make the energy gap only polynomially small 
A cx e/A 8 23 if an boost and projection Hamiltonians 
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FIG. 1: The same as Fig. (2) in this figure shows, in a 
complicated circuit, how the CNOT interacting qubits is mod- 
ified by inserting teleportation boxes on each qubit's upstream 
and downstream part, so that the energy gap is only polyno- 
mially small. Each dot represents a row of two quantum dots, 
label I stands for identical transformation Hamiltonian h(I), 
H for Hadamard transformation Hamiltonian h(H), and P(0) 
for projection Hamiltonian h(\0),X). 

have the same amplifying factor A. To determine mag- 
nitude of A, one only needs to count the total number 
of qubits in the circuit, say L, which is proportional to 
the number of control operation in an algorithm, then 
the probability of finding all electrons on final rows is 
P « (1 — C/X 2 ) L with C being 8, the maximum length 
of qubit. In order to have appreciable P, we set A ~ L 1 / 2 , 
hence A oc e/L 4 . The details can be found in [T3 |. 

C. Energy Gap When Projecting Small Fraction of 
a State 

In the previous section the rule for finding scale of the 
energy gap is under the assumption that when a projec- 
tion Hamiltonian h(\~/),X) is applied, H/\/H 2 + \b\ 2 is 
appreciable for the ground state on row just before the 
projection Hamiltonian: 

o| 7 >+6|7>, (8) 

where I7) = |0) (|1)) and (7) = |1) (|0». The ground- 
state wavefunction concentrates on the last row, hence 
the first excited state wavefunction cannot have appre- 
ciable weight there because otherwise (ipi |^o) 7^ 0. When 
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FIG. 2: A six-row single qubit ended with the projection 
Hamiltonian h(\0),X). 
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FIG. 3: Energy gap A verse A with h(\0),X) applied on the 
last row of a 6-row single qubit, where h = 10e(/ + aa z — 
Vl — a 2 a x ). From top to bottom, lines correspond to a = 
0, 0.9, 0.99, 0.999, 0.9999, 0.99999. 



evaluate 1/x on a qubit, the projection Hamiltonian on 
the qubit itself doesn't contribute to For example, 
concerning a single qubit, as shown in Fig.®, with only 
identical transformations h(I) and ended by h(\Q),X), if 
h° = E(I — a x ) so that |^o) on the first row is |0) + |1), 
then the energy gap A is almost independent of A, as 
shown in the top line of Fig. © . 

However, if in Eq.@ M/VM 2 + W < 1, then A 
depends on A until A reaching \/\a\ 2 + |6| 2 /|a|. This 
is because when A < ^/|a| 2 + |6| 2 /|a|, the ground-state 
wavefunction has little weight on the last row, and the 
first excited state concentrates there, hence 1 jx is small, 
leading to small energy gap. When A > \J\a\ 2 + |6| 2 /|a|, 
ground state wavefunction has large part on the last row, 
then just like the above situation, energy gap is not fur- 
ther affected by increasing A. 

To confirm the above analysis, we numerically calcu- 
late the energy gap of a 6-row single qubit ended with 
projection Hamiltonian, as shown in Fig- (EJ - The bound- 
ary Hamiltonian is h° = 10e(7 + aa z — \/l — a 2 a x ), 
all other Hamiltonians except for that at final row are 
h? (I) with j = 1, 2, 3, 4, and on the final row 
there is a projection Hamiltonian h(\0),X). By tun- 
ing a, we can determine what fraction of wavefunc- 
tion is projected from the 5th row to the last row. At 
a = Q, 0.9, 0.99, 0.999, 0.9999, 0.99999, on the 5th row 
the ground state wavefunctions are a|0) + 6|l) with a/6 = 
1, 0.23, 0.071, 0.022, 0.0071, 0.0022. Fig.® shows that 
the energy gap is A oc e/A 2 as A < |-\/|a| 2 + |6| 2 /a|, and 
when A > |-\/|a| 2 + |6| 2 /a|, A becomes independent on A. 
The independent A is proportional to e|a/i/|a| 2 + |6| 2 | 2 . 
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In order to make the ground-state wavefunction con- 
centrate on the last row so that measurement corre- 
sponds to the desired state, A must be larger than 
IVM 2 + \b\ 2 /a\, Thus the energy gap is determined by 
the fraction of state been projected. If \a/ \/\a\ 2 + \b\ 2 \ is 
exponentially small, which may happen in certain case, 
then the energy gap is exponentially small. Fortunately, 
this doesn't happen to the GSQC implement of Quantum 
Fourier Transform, there all projection Hamiltonians are 
applied to teleportation circuit, and \a/b\ — 1. However, 
it plays a role in the algorithm presented in the following 
section. 

For multiple interacting qubits, if \a/ \/\a\ 2 + \b\ 2 \ <C 1 
in Eq.JSJ), the rule of finding energy gap needs modifica- 
tion: With all qubits ended with either a projection or a 
boost Hamiltonian containing the same amplifying factor 
A ^> 1, when estimating 1/x for any qubit, say qubit A, 
(i) at first x is set to 1; (ii) the boost Hamiltonian, or 
the projection Hamiltonian, on qubit A itself increases x 
by multiplication of A or Min(A, |i/|a| 2 + |6| 2 /a|); (Hi) if 
qubit A directly interacts with another qubit, say qubit 
B by Hamiltonian ft. as, then we determine, excluding 
qubit A, on the qubit B the amplitude ratio of the up- 
stream part (divided by Hab) over its final row, 1/xb, 
and contribution to 1/ xb from other qubits are found one 
by one according to Eq. © : if the upstream part of qubit 
B doesn't coexist with the states on final rows of a qubit, 
except for qubit A, then xb should be multiplied by A 
(ended with boost Hamiltonian) or X\a" /W\a"\ 2 + |6"| 2 | 
(ended with projection Hamiltonian); (iv) finally, the 
value of 1/x on qubit A should be multiplied by 1/xb 
or Hil/x B if more than one qubit directly interact with 
qubit A. 

It is easy to find that when \b/a\, \b"/a"\ w 1 and 
A > 1, we get the same result as the previous subsection. 
After 1/x's on all qubits being evaluated, the minimum 
1 jx gives the energy gap scale as 

A cx e(l/x) 2 min . 



III. QUANTUM ALGORITHM BY GSQC 

There are some interesting properties for the GSQC. 
Although it was shown^^ that, concerning on time cost, 
a quantum computer composed of (time varying) local 
Hamiltonians is equivalent to standard circuit quantum 
computer, GSQC provides some insights to design quan- 
tum algorithm for certain problems. For example, the 
projection Hamiltonian, which corresponds to measure- 
ment in standard paradigm, can amplify the probability 
at a particular state. Here we are not claiming that the 
GSQC is more powerful than standard quantum com- 
puter, however, the GSQC does provide a direct approach 
for certain problem, as shown below is the algorithm for 
the SAT problems. 

At first we give the simplest example, considering that 
qubit i CNOT controls an ancilla qubit that is at the 
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FIG. 4: A filter for the clause i = 0. 



pin 



I0> 






1 4 


1 • 




I 

» 4 

CNOT ( 


► 4 

i j 

> 4 

CNOT 




> 4 




B 

1 4 


B 

i i 



ancilla qubit i J 



FIG. 5: A filter for the clause i + j — 1. 

right side in Fig.Q, and their boundary Hamiltonians 
make the ground state on their first rows are |0) + |1) 
and |0), respectively. On last rows the ground state is 
|0)|0) + If we apply a boost Hamiltonian on qubit 

i and a projection Hamiltonian ft(|0),A) on the ancilla 
qubit, then at the ground state the state on final rows 
becomes |0)|0). The large value of A makes sure that 
there is large probability to find two electrons on the 
final rows of the two qubits at the ground state. So by 
choosing projected state on the ancilla qubit, we can have 
the selected state |0) on qubit i, and prevent the other 
state |1) from reaching its final row. If qubit i entangles 
with other qubit, such as |0) |a) + 1 1) \f3) , the entanglement 
of |0)|a) will not be affected. Thus we call circuit in 
Fig.(0} a filter for the clause i = 0. 

Another example makes more sense. Lets consider a 
SAT problem with clauses, each of which involves two 
qubits, say qubit i and j, and requires i + j = 1. We can 
implement this clause by the GSQC circuit in Fig.j5j|. In 
this figure there are three qubits: qubit i, qubit j and an 
ancilla qubit that is at the left side in the figure. It's easy 
to find that if on the first row \i) — cti\0) + /3i\l), \j) = 
ctj\0) + (3j\l) and the ancilla qubits at |0), then at the 
ground state on the final rows of the three qubits the state 
is \i)\j)\ancilla) = (fta J -|l)|0)+Q!^|0)|l))|l) ) which sat- 
isfies the clause. Thus circuit in Fig.JSJ filters out states 
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FIG. 6: A GSQC circuit solving SAT problem with each 
clause involving several bits. Box labeled "Filter Box" repre- 
sents filter like Fig. gj (one-bit clause), Fig. JSJ (two-bit clause) 
or Fig. (three-bit clause). 

not satisfying this simple clause and lets through those 
satisfying states. It is important to note that at the be- 
ginning if the satisfying states entangle with other qubits 
not showing in the figure, these entanglements keep un- 
touched. 

The property of GSQC brings up new quantum algo- 
rithm naturally. Here we present one to solve the SAT 
problems as shown in Fig. © , a GSQC circuit to solve a 3- 
SAT problem with only 9 bits. It's easy to be extended to 
iV-bit iC-SAT problems. Each clause is implemented by 
a "filter box" , and the circuit inside each filter box makes 
sure that on rows immediately below it the ground state 
satisfies the clause Ci, or we can say those unsatisfying 
states are filtered out. This can be realized by projection 
and boost Hamiltonians like in Fig.Q) and Fig.©. 

In Fig.©, the initial state on the top rows of qubit 
from 1 to 9 is (|0) + |1»(|0) + |1))...(|Q) + |1)), which 
is enforced by the boundary Hamiltonians, h = E(I — 
<j x )\ the clause involving qubit 1, 2 and 3 is implemented 
by filter box 1, the clause involving qubit 2, 3 and 4 
implemented by filter box 4, the clause involving qubit 
3, 4 and 8 implemented by filter box 6, etc. 

When all constraints are implemented, at ground state 
the states measured on the final rows of the N qubits 
should be superposition of all states satisfying all con- 
straints. No backtracking is needed. 

Concerning energy gap, unlike the circuit for quantum 
Fourier transform, in which the energy gap is determined 
by the number of control operation [lj], the SAT prob- 
lems is more complicated to evaluate because it might 
involve the situation to project a very small fraction of 
state as shown in section III CI For example, if one con- 
structs a GSQC for the Grover's search problem with one 
condition to find a unique satisfying state from 2 N states, 
then he will find that there is an ancilla qubit containing 
such unnormalized state 

2 N -1 

|0) (satisfying) + |l)|unsatisfying (i) ) (9) 

i=l 

before the projection Hamiltonian h(\0), A). In order to 



amplify the amplitude of the correct state on the final 
row, it requires A > 2 N ' 2 . Its energy gap is hence less 
than 2~ N , which is consistent with the limit set by many 
other works H [H [13. 

IV. EXAMPLE: THE 3-BIT EXACT COVER 
PROBLEM 

Up to now the filters, Fig.(gJ and Fig.©, we have 
given are trivial, and now we give an example on how 
to implement a filter for a serious problem. We focus 
on the 3-bit Exact Cover problem 6] , an instance of SAT 
problem, which belongs to NP-complete. Following is the 
definition of the 3-bit Exact Cover problem: 

There are N bits zi, Z2, zn, each taking the value 
or 1. With O(N) clauses applied to them, each clause is 
a constraint involving three bits: one bit has value 1 while 
the other two have value 0. The task is to determine the 
N-bit state satisfying all the clauses. 

A. GSQC Circuit for the 3-bit Exact Cover 
Problem 

The algorithm is implemented by the circuit in Fig.©. 
Each filter box, in our algorithm, involves three qubits, 
say qubit i, j and k, which are represented by gray dot 
columns in Fig.Q. We add two ancilla qubits: qubit 1 
and qubit 2, which are represented by dark dot columns. 
Qubit i, j and k at the first row are in the state (|1) + 10)) 
if they have not experienced any clause yet, and the two 
ancilla qubits are in the states |0) and |0) on top rows 
by selecting proper boundary Hamiltonians, where I7) 
corresponds to the state of ancilla qubit 1, and (7) to the 
state of ancilla qubit 2. 

Inside the dashed triangle of Fig.J7J), after the first 
CNOT, we obtain state + |6)|0); after the second 

CNOT: |1)|1)|0) + |6)|0)|0) + |6)|1)|1) + |i)|0}|l); after 
the third CNOT: 

|i)(|i)|o)|o) + |o)|i)|o) + |o)|o>|i) + |i)|i)|i)) 
+ |6)(|i)|i)|o) + |o)|i>|i) + |i)|o)|i) + |o)|o)|o)). 

Immediately below the triangle, if the system stays at the 
ground state, if electron in ancilla qubit 1 is measured 
to be on the row labeled by X and at state |1), and if 
the three electrons on qubit i, j, k are all found on the 
rows labeled by X, then the three-qubit state satisfies 
the clause except for |1)|1)|1). 

The ancilla qubit 2, starting at state |0), experiences 
CNOT gates controlled by qubits j and k, and R(±tt/4) 
transformations, defined in 18] as i? !/ (±7r/4), as shown 
within the dotted pentagon in Fig.J7J. All those transfor- 
mations happened inside the dotted pentagon are equiv- 
alent to a Toffoli gate except for some unimportant 
phases 18]: if both qubits j and k are in state |1), then 
the ancilla qubit 2 reverses to state |1), otherwise, it 
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FIG. 7: A filter for the clause i + j + k = 1. The labels on 
the lines stand for corresponding Hamiltonians: / for 
CWOT for h(CNOT), P(l) for projection A) et. al. At 

the final rows, B/I represents boost Hamiltonian h(B, A) if 
there is no more clause to be applied to this qubit, otherwise, 
represents identical transformation Hamiltonian h(I). There 
are teleportation boxes, not shown in figure, inserted on all 
qubits between two control Hamiltonians. Some dots marked 
by X or Y are for demonstration convenience in text. 



remains at state |0). After this nearly Toffoli transfor- 
mation, if at ground state electrons in qubit j, k and 
ancilla qubit 2 are found on rows labeled by Y, and if 
ancilla qubit 2 is at |0), then the three qubits will be at 
|0)(|0)|0) + |1)|0) + |0)|1)). Thus if at ground state all 
electrons are found on rows immediately below both the 
dashed triangle and the dotted pentagon, and if ancilla 
qubit 1 is at |1) and ancilla qubit 2 at |0), then the three 
qubits i, j, k satisfy the clause: 



|i)|5)(|i)|o)|o) + |o)|i)|o) + |o)|o)|i)) 



(10) 



In order to make the satisfying states pass through 
the filter box with large probability, we add projection 
Hamiltonians and boost Hamiltonians as shown in the 
lower part of Fig.Q. The projection Hamiltonians on 
final rows of the two ancilla qubits limit and amplify the 
amplitude of the states we prefer: ancilla qubit 1 at |1), 
and ancilla qubit 2 at |0). If a qubit does not experience 
any more clause, it will end with a boost Hamiltonian, 



otherwise, its quantum state will be teleported to a new 
qubit through teleportation box, not shown in Fig.lJZJ), 
and the new qubit experiences more clauses. Thus the 
projection Hamiltonians on two ancilla qubits and boost 
Hamiltonians on the three qubits make sure that the 
ground-state wavefunction concentrates on the final rows 
in Fig. 10 with state at Ea. (flT))l . 

Noting that in the filter box all the three qubits i, j, 
and k always act as control qubits, thus the entangle- 
ments of these three qubits with other qubits not involved 
in this particular clause still keep the same. When adding 
a clause, the resulted states satisfying this clause will also 
satisfy all previous applied clauses. Thus unlike classical 
algorithm, no backtracking is needed. 



Energy Gap Without Projecting Small Fraction 
of State 



In this subsection, we assume applying each clause 
does decrease the number of satisfying state gradually, 
or equivalently, the projection Hamiltonian in the two 
ancilla qubits in each filter box, Fig-.Q, does project 
appreciable part of state on the second last row. This 
assumption may not be correct in many SAT problems, 
especially close to a c . 

In the circuit of Fig.JHJ, if there is at least one solution, 
and all electrons are simultaneously found on the final 
rows of all qubits, then the reading of the iV-bit state 
satisfies all clauses. 

In order to keep the energy gap from being too small, 
like in |14| . on every qubit teleportation boxes are in- 
serted between two control Hamiltonians, thus the to- 
tal number of qubits increases while the energy gap 
A oc e/A 8 if all the boost and the projection Hamilto- 
nians have the same value of amplifying factor A. 

For one clause, or a filter box, it needs 10 teleporta- 
tion boxes (each teleportation box adds two more qubits) 
on the original five-qubit circuit, noting that on the end 
of qubit i, j and k in Fig.Q teleportation boxes are 
needed because more clause will be added. Thus adding 
one more filter box means adding 20 more qubits. The 
number of clause for a NP hard 3-bit Exact Cover prob- 
lem is about the same order as the number of bits iV0, 
say aN with a being O(l), then there are about 20aN 
qubits and each of them ends with either a projection or 
a boost Hamiltonian. Probability of finding all electrons 
at the final rows is approximately 



p m (i ~c/\ 2 y 



(ii) 



where C = 8, the length of the longest qubit 14] . It is 
assumed that, at ground state, in each filter box the an- 
cilla qubit 1 and 2 have appreciable probability in 1 1) and 
|0) states, respectively, before projection Hamiltonians. 
Later we will address the situation when this assumption 
is violated. 

In order to make the probability independent of num- 
ber of bits N, we take A 2 = DN, where D is an arbitrary 
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number. Then as N becomes large, we obtain 



P w (1 - C/{DN)) 20aN « e~ 
and energy gap is 01 

A oc e/A 8 oc e/(D 4 N 4 ) 



20aC/D 



(12) 



(13) 



from which one can estimate time cost. 

To make the GSQC circuit at ground state, we can 
use adiabatic approach: first we set A = 1 for boost and 
projection Hamiltonian on final rows of all qubits, and 
replace the single qubit Hamiltonian between the first 
two rows of all qubits by a boost Hamiltonian 



ti(B,\') = e 



^C\d + CtC 2 - y, (c\C2 + h 



A' 



GL4) 



so that the wavefunction amplitude of the first row is 
boosted as A' S> 1. Now in the ground state the elec- 
trons concentrate at the first rows as 1/A' — > 0, thus 
the ground state is easy to be prepared, and the energy 
gap A oc e/n 2 with n = 8 being the length of the longest 
qubit. The next step is turning the quantity 1 /A' to 1 adi- 
abatically, during which the energy gap remains at e/n 2 
and the ground-state wavefunction spreads to other rows 
from the first row. The third step is turning 1/A from 1 
to l/y/DN adiabatically. In this process the energy gap 
decreases monotonically from e/n 2 to what we obtained 
above: e/D 4 N 4 , and the ground-state wavefunction con- 
centrates on the final rows of all qubit as we wish. Thus 
the scale of time cost is about T oc 1/A 2 oc A 8 0, local 
adiabatic approach may reduce the time cost further pflj. 



C. Energy Gap for SAT Problems 

Above analysis is under the assumption that the num- 
ber of satisfying states gradually decreases as the clauses 
are implemented one by one. There is a situation that 
might hurt our algorithm: after adding one more clause, 
the number of satisfying states drops dramatically. Just 
like what happens to Grover's search algorithm, in which 
the number of satisfying states drops from 2 N to 1, and 
as shown in Eq.@, our algorithm involves a projection 
Hamiltonian on an ancilla qubit to project an exponen- 
tially small fraction of a state, thus the energy gap eval- 
uation in the above subsection becomes invalid. 

Does this happen to the general SAT problems? In 
it was suggested that close to the threshold a c com- 
putational complexity might be related with the form- 
ing of a backbone, each of a subset of bits has average 
value close to 1 or in the subspace of satisfying states. 
The existence of the backbone means that most satisfy- 
ing states contain the state represented by the backbone, 
and if adding one more clause kicks out the states con- 
sistent with the backbone from satisfying subspace, the 
number of satisfying states drops dramatically, and this 
corresponds to projecting a small fraction of state. 



Performance of our algorithm is not affected by form- 
ing of backbone, however, as more clauses applied, the 
disappearance of the already existed backbone in the sat- 
isfying subspace surely hurts. There is a criterion deter- 
mining efficiency of our algorithm: the ratio Sj/Sj+i, 
with Sj being the number of solutions when the jth 
clause is applied, and Sj+i the number of solutions when 
the (j + l)th clause is applied. For example, So/ Si = 8/3 
for 3-bit Exact Cover problem. If Sj/Sj+i ^> 1, on 
the ancilla qubit of the (j + l)th filter box, the prob- 
ability of finding electron on its final row will be p rs 
(1 — CSj/(\ 2 Sj+i)). To make sure of appreciable proba- 
bility of finding all electrons on the final row of all qubits, 
an overhead factor J Sj/Sj+i for A on the ancilla qubit 
is needed, hence the amplifying factor in the projection 
Hamiltonian on the ancilla qubit should be \yj Sj/Sj+±. 
According to the analysis in Sec lII CI the energy gap 
might be also determined by the parameter Sj/Sj+i. Be- 
cause in a filter box, the ancilla qubit will end after the 
projection Hamiltonian, which should be at the position 
of qubit 8 or qubit 10 in Fig.JTJ without the dotted line 
following. According to the rule described in section lTl CI 
the parameter 1 / x on this ancilla qubit should be 



(15) 



\ 2 Min A 



The energy gap thus is 



Min I 



eS 



A 4 S* 7 - 



(16) 



If this ratio Sj+i/Sj happens to be exponentially small, 
then our algorithm cannot solve the SAT problem in 
polynomial time. We cannot know in advance what 
Sj+i/Sj is, however, we might be able to identify back- 
bone by trials, and then choose proper order to imple- 
ment clauses so that Sj+i/Sj always can be kept not 
too small. However, if the NP-Complete problem means 
that one can never avoid an exponentially small Sj+i/Sj, 
then the quantum algorithm cannot solve NP-Complete 
problem in polynomial time. 



V. CONCLUSION 

In conclusion, we have demonstrated that a ground 
state quantum computer can solve a general SAT prob- 
lem. A specific example, the 3-bit Exact Cover problem, 
is given. We show that a 3-bit Exact Cover problem can 
be solved by the quantum algorithm described here, and 
the time cost is related with the number of bits N and the 
parameter Sj+i/Sj. If Sj+i/Sj stays only polynomially 
small, then the presented algorithm can solve this SAT 
problem in polynomial time. It will be interesting if one 
finds the equivalent algorithm by standard paradigm. 
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